Retarded Casimir-Polder force on an atom near reflecting microstructures 
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We derive the fully retarded energy shift of a neutral atom in two different geometries useful for 
modelling etched microstructures. First we calculate the energy shift due to a reflecting cylindrical 
wire, and then we work out the energy shift due to a semi-infinite reflecting half-plane. We analyze 
the results for the wire in various limits of the wire radius and the distance of the atom from the 
wire, and obtain simple asymptotic expressions useful for estimates. For the half-plane we find an 
exact representation of the Casimir-Polder interaction in terms of a single, fast converging integral, 
which is easy to evaluate numerically. 

PACS numbers: 31.70.-f, 41.20.Cv, 42.50.Pq 



I. INTRODUCTION 

The explosive rate of developments in nanotechnology 
as well as in the manipulation of cold atoms has meant 
that interest in atom-surface interactions has increased 
strongly in recent years. What were once tiny, elusive ef- 
fects are now dominant interactions, or, as the case may 
be, a major nuisance in some experimental set-ups. Mo- 
tivated by a common type of microstructure, which con- 
sists of a protruding ledge fabricated by successive etch- 
ing and possibly a thin electroplated top layer, we have 
recently studied the force on a neutral atom in close prox- 
imity of reflecting surfaces of either cylindrical geometry 
or that of a semi- infinite half-plane |Jj] . In the absence of 
free charges or thermal excitations, the interaction of the 
atom with the microstructure is dominated by Casimir- 
Polder forces , which are due to the interaction of the 
atomic dipole with polarization fluctuations excited by 
vacuum fluctuations of the electromagnetic field. If the 
atom is sufficiently close to the surface of the microstruc- 
ture, then the interaction between the atomic dipole and 
the surface is purely electrostatic and retardation can be 
neglected, which was the case investigated in Ref. 
Then one does not need to quantize the electromagnetic 
field, but can work with the classical Green's function of 
Poisson's equation. The only difficulty lies then in the 
geometry of the problem. 

However, in experimental situations one more often 
finds that retardation is in fact important, as the dis- 
tance of the atom from the surface of the microstructure 
is often commensurate or larger than the wavelength of 
a typical atomic transition. This is the case we investi- 
gate here, again for microstructures of two types of ge- 
ometries: a cylindrical reflector of radius R and infinite 
length, and a reflecting half-plane. 

Various versions of this problem have been studied be- 
fore, both analytically and numerically. Probably the 
first to consider the interaction between an atom and a 
metallic wire, according to was almost 75 years ago 
Zel'dovich j^. This problem was then revisited and ex- 
tended by Nabutovskii et al. [5|], and subsequently by 
Marvin et al. [6l|. In Nabutovskii's paper a dielectric 



cylinder is envisaged to be surrounded by a cylindrical 
shell of vacuum which in turn is surrounded by a rar- 
efied gas of polarizable particles. The interaction energy 
of a single particle is then calculated through the work 
done by the force (obtained from the stress tensor) due 
to the fluctuating electromagnetic fields, in the limit of 
zero density of the surrounding gas. The asymptotic re- 
sults obtained there (Eq. (23) and Eq. (24) of Ref. [| ) 
are, according to Ref. 3], valid only for dilute dielectric 
materials; they diverge in the perfect-reflector limit. 

On the other hand, the work by Marvin et al.Q, mo- 
tivated by ^, ^ and based on a normal-mode expansion 
and a linear-response formalism Q, gives the same gen- 
eral formula for the interaction between a point parti- 
cle and a cylinder [their Eq. (4.10)] as the equivalent 
result in fE]. We have no reason to believe that the re- 
sult in [6] is incorrect in the perfect-conductor limit, as 
it reduces to our previous result [l| in the electrostatic 
limit. Moreover, Ref. [6 manages to recover the original 
Casimir-Polder result [2| in the large-radius limit of the 
cylinder. This suggests that the general expression in Q 
is probably correct, only that the perfect-conductor limit 
does not commute with the asymptotic limit of the zero 
radius (or large distance of the atom from the cylinder) 
studied there. In the small-radius limit, the result for the 
interaction between an atom and a metallic filament, in 
both retarded and non-retarded limits, is also given by 
ffl- , 

Marvin et al.'s work [6[ is certainly the most compre- 
hensive, but due to its generality it is also quite cum- 
bersome to apply, which is mainly done numerically for 
just a few examples [13] • Further numerical studies of 
the interaction of atoms with macroscopic cylinders can 
be found in Refs. dH, [H, [H, [11 . 

By contrast, in this paper we are after a relatively sim- 
ple theory that allows one to estimate the force between 
an atom and a cylindrical reflector at any distance and 
cylinder radius. To this end we are not interested in the 
precise dependence of the interaction on material con- 
stants of the reflector, and therefore we work with the 
model of a perfectly reflecting surface. 

As discussed in Ref. [l| , we also determine the force be- 
tween an atom and a semi-infinite half-plane, in order to 
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facilitate estimates for common types of microstructures 
that consist of a ledge protruding from a substrate. The 
Casimir-Polder interaction between an atom and such a 
half-plane has also studied before, but only in the ex- 
treme retarded limit of very large distances of the atom 
from the surface [HI . To the best of our knowledge no for- 
mula for the interaction in the intermediate region, when 
the distance of the atom from the surface is comparable 
to the typical wavelength of an internal transition in the 
atom, has been derived yet. Recent work of Mendes et 
al. [la], dealing with wedges, does not include the gen- 
eral result in the half-plane geometry as a limiting case 
of a zero-angle wedge. 



dipole approximation, which simplifies the expression for 
the energy shift to 



2£o 



Eh + UJX 



(6) 



where we have introduced the abbreviation Eji = Ej—Ei. 
For brevity and presentational clarity we shall henceforth 
also abbreviate the matrix elements of the atomic dipole 
moment as 



\|^\^\{J\^^\^)\ 



(7) 



II. FIELD QUANTIZATION AND THE 
ENERGY SHIFT 



III. ENERGY SHIFT NEAR A PERFECTLY 
REFLECTING WIRE 



The complete system of an atom interacting with 
the quantized electromagnetic field is described by the 
Hamiltonian 



H = FAt, 



H-p; 



Field 



Int 



(1) 



We choose to work with /i-E coupling, i.e. our interaction 
Hamiltonian is 



Hint — H ■ E 



(2) 



Quantization of the electromagnetic field is done by way 
of a normal-mode expansion of the vector potential in 
terms of photon annihilation and creation operators for 
each mode A and polarization a, 
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A,CT 



(3) 



To describe a mode we use the composite index A instead 
of a wave vector, as we shall be working in cylindrical co- 
ordinates where the quantum number of the azimuthal 
part of the mode function is discrete, but the other two 
are continuous. We work in Coulomb gauge, V • A(r), so 
that the normal modes F(r) satisfy the Helmholtz equa- 
tion. 



iF(r) . 



(4) 



The energy level shift due to the interaction ^ can be 
calculated perturbatively. For our system in state |«;0), 
i.e. the atom in state \i) and the electromagnetic field in 
its vacuum state |0), the lowest non- vanishing order of 
perturbation theory is the second, so that 



AW = 



A 



-/i-E|i;0 



Ei - {Ej + wa) 



(5) 



As the relevant field modes can be expected to vary 
slowly over the size of the atom, we make the electric 




FIG. 1: Atomic electric dipole moment in the vicinity of a 
perfectly reflecting cylinder of radius R. The normal modes 
F^'^'(x) in this geometry are given by Eqs. (|13p and (|14p . 



First we wish to calculate the energy shift of an atom 
near a perfectly reflecting and infinitely long cylindrical 
wire of radius R. It is advantageous to work in cylindrical 
coordinates, cf. Fig. [T] 

In order to find two independent transverse vector field 
solutions of Eq. ([I]), we make use of the re pre sentation 
theorem for the vector Helmholtz equation 10.411]. 
If $(x) is a solution of the scalar Helmholtz equation then 
the two independent solutions of the vector equation are 
given by 

F(i)(r)-(Vxe,)ci>(r) , (8) 
F(2)(r) ^ iy X (V X e^) $(r) . (9) 

The particular choice of the constant unit vector is 
motivated by the symmetry of our problem and lets us 
to identify the solutions F(^)(r) and F(^^(r) with the 
transverse electric (TE) and transverse magnetic (TM) 
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modes, respectively. In cylindrical coordinates the scalar 
Helmholtz equation has the solutions of the form 

<P,z) = N[ cos6„Mkp) + smS,nY^{kp) ] e'^^+i- 

(10) 

where Jm{kp) and Y^(kp) are Bessel functions of the first 
and second kind [la §9]. The separation constants sat- 
isfy = fc^ + K^, and m is an integer. The phase shifts 
5m describe the superposition of regular and irregular 
solutions. In free space only regular solutions Jm{kp) 
are admissible, and 6m = 0. In the presence of the per- 
fectly reflecting wire, the phase shifts serve to make the 
electromagnetic fields satisfy the boundary conditions on 
the surface of the wire. The normalization constant N is 
chosen such that 



d3rFiT^*(r).Fi-)(r)=5„ 



d{k-k') 



(11) 



is met. Setting cos 5m = 1, s'm 5m = 0, one can derive 
quite easily that N — {2Trk)~-^. 

On the surface of a perfect conductor, the tangential 
components of the electric field and the normal com- 
ponent of the magnetic field vanish. Therefore, at the 
surface p — R of the cylindrical wire we must have 



= E, and B„ 



0. These boundary conditions 



determine the phase shifts as 



tan^™ = 



YLikR) 



, tan<5™^ 



JmjkR) 
YmikR) 



(12) 



and (IIOII , the normalized mode 



^"^'(r), A = {fc,TO, k}, that satisfy the bound- 
R are given by 



According to Eqs. 

functions F 
ary conditions at p 

"im Jmikp)Y^,ikR) - Ym{kp)J'm{kR) 



1 
2^ 



kp 



^/iKm+Y^^) 



(13) 



j'^{kp)Y:^{kR) 



Y:,,{kp)j'^{kR) „ 



VJiikW+Klm 



{P,(l},z) 



1 

2^ 



IK J'^{kp)Ym{kR) - Y;^{kp)Jm{kR) 



ruK Jm{kp)Ym{kR) - Ym(kp) JmjkR) 
'^P ^Jl{kR)+Y^l{m 
fc Jm{kp)Ym{kR) - Ym{kp)Jm{kR) „ 
^JUkR) + Y^{kR 



(14) 



These mode functions can now be substituted into Eq. ([6|) 
for obtaining the energy shift of an atom positioned at 
r = (p, 0, z). However, what we want to calculate here 
is only the correction to the energy shift caused by the 
presence of a perfectly conducting surface, rather than 
the whole energy shift due to the coupling of the atom 
to the fluctuating vacuum field, which would include the 
free-space Lamb shift. Therefore we need to subtract 
the energy shift caused by the vacuum fluctuations of 
the electromagnetic field in free space, which is obtained 
by either letting the phase shifts 5m — > or equivalently 
taking the limit i? — s- 0. In the limit of vanishing radius R 
of the cylinder the behaviour of the mode functions (|13p . 

is dominated by the singular behaviour of Ym{kR) 
and Y^^(kR) at the origin, which causes the phase shifts 
p2)) to vanish. The rcnormalized energy shift Aiy''" = 
AW — lini/j^o is found to be of the form 



47r£o ^ 



with 



(15) 



■E' 

m— 



„ ^ \ ' / 
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m \ ' 

Vp) 



{Jm{kp)Y:^{kR) - Ym{kp)JUkR)y 
J'i{kR)+Y:^{kR) 



Jrnikp) 



{j'^{kp)Ym{kR)-Y;^{kp)Jm{kR)) 



Jl{kR)+Y.l{kR) 

{j:^{kp)Y:^(kR) - Y:^{kp)j^{kR)Y 

J2{kR)+Y:;}{kR) 
' {Jm(kp)Ym{kR) - Ym{kp)Jm{kR)) 



- &p) 



Jl(kR)+Y^{kR) 

fcy r {Jmikp)YmikR) ~ Ym{kp)Jm{kR)f 
Lo) Jl{kR) + Yl{kR) 



Jmikp) 



Jlikp) 



(16) 



(17) 
(18) 



4 



where the primes on the sums indicate that the m = 
term is weighted by an additional factor of 1/2. It ap- 
pears that the k integrals fail to converge, but this is a 
common feature in such calculations caused by the dipole 
approximation, see e.g. jj^]. As we shall see, convergence 
is in fact brought about by the Bessel functions, which 
come to bear if the k integral is replaced by an integral 
over to = y/ + k?. 

As the Bessel functions Jm{x) and Ym{x) are both 
oscillatory for large x, we wish to rotate the integra- 



tion contour in the complex k plane, in order to get an 
integrand that is exponentially damped for large argu- 
ments. To this end we introduce the Hankel functions 
H'^k\x) = J^{x)+iY^{x) and hI^\x) = [H'^k\x)]* = 
Jm{x) — iYm{x), in terms of which wc can rewrite the 
energy level shift in such a form that there are no poles 
in the first quadrant of the complex k plane, as is re- 
quired for the rotation of the integration contour. This 
step greatly simplifies further analysis. 



= -Re - E ' Tdfc k Hdn {-,[H'i'\kp)f^^^^^ + ^[H^^Hkp)?^^j^] , (19) 

E,^-Re^f2' Hdk k Hdn I [H^W ^kp)? + J!^4[H(i)(fc,)]2 1 , (20) 

^ ^0 Jo Jo Ej,+ij[ H'A \kR) k^P^u;^ H}n'{kR)} 

= -Re 2 g ' r dfc k r dn — ^ 1 4 [H^? (Ml^ I ' (^l) 



We now transform the k integration in Eqs. (fT9|) -([2T |) 
into an integration over uj — \J -|- A:^, and note that 
on the interval < < k the integrands become pure 
imaginary and therefore do not contribute if added to the 
real part of the integral. We can therefore shift the lower 
limit down to the origin 

/"OO /"OO 

/ dw — > / dw (22) 

J K Jo 

without affecting the result. Further, we note that the 
functions (z) and H.'m'' [z) have no zeros in the first 



quadrant of the complex plane 18, Fig. 9.6], so that the 
contour of the w-integration can be rotated from the pos- 
itive real to the positive imaginary axis, to — s- \lo. Then 
the oscillatory Bessel functions turn into the modified 
Bessel functions according to 18, 9.6.3 & 5] 

J™(iz) = e'™^/2/™(z), (23) 
ilW(iz) = --e— /2if™(z). (24) 

TT 

Taking the real part and going to polar coordinates, 
where the angle integrals are elementary, we find that 




Note that the effect of our manipulations has been that the integration variable k in Eqs. (P5)) - (|?7|) has been ro- 
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tated by 7r/2 in the complex plane compared to Eqs. (|19p- 

The final result for the energy shift, Eq. (fT5|) with 
Eqs. (P5)) - (P7)) . is a sum over a series of rapidly converg- 
ing integrals, which, unhke Eqs. P^ - (fT5)) . is reasonably 
easily evaluated numerically. However, as the functions 
'^p.(t>.z{Eji,d,R) are quite cumbersome and it is not pos- 
sible to find exact closed form expressions for them, we 
now look at their asymptotics in various limiting cases, 
which is very useful for analytical estimates. 



A. Asymptotic regimes 

There are three length scales in the problem: the 
distance of the atom from the surface of the cylin- 
der d — p — R, the radius of the cylindrical wire R, 
and the wavelength of a typical transition in the atom 
\ji oc i/Eji. Accordingly we get six different asymptotic 
regimes, three non-retarded and three retarded. The cri- 
terion as to whether retardation matters is the relative 
size of the distance d of the atom from the surface and the 
wavelength Xji of a typical transition: if the atom is very 
close to the surface then its interaction with the surface 
is entirely electrostatic [1] , whereas retardation begins to 
play a role once d ^ Xji or larger, because then the in- 
ternal state of the atom is then subject to non-negligible 
evolution during the time a virtual photon mediating the 
interaction would take to travel from the atom to the 
surface and back. First we shall deal with the three non- 
retarded cases, and then with the three retarded ones. 



electrostatic limit {Eji 0) and the large-radius limit 
(i? — > oo) of the energy shift commute. 

The limit of large radius was studied in great detail 
in . Application of the summation formula derived in 
Appendix A of [g'] to Eqs. ([251)- (H?]) leads to the original 
Casimir-Polder result [2] for the interaction between an 
atom and a plane, perfectly refiecting mirror: 



1 



27rd3 
1 

2^ 



dr] 



drj 



g-2d_Bjir) 
(1+772)2 ' 

(1 + 772)2 1+772 



(29) 
(30) 



If we now take Xji to be much greater than d, we repro- 
duce the result (p8)) of the previous section. 



3. i? < d < A, 



In this case we again start by taking the limit Eji 
in Eqs. (|25| - ([27|) and obtain the electrostatic expression 
derived in [l| . In the limit of the radius of the wire be- 
ing much smaller than the distance d, the energy shift is 
dominated by summand with lowest m in Eqs. (|25 p ~(|27 p 
[ij. Asymptotically one gets 



1. d<i?<Aji 

If Xji is larger than any other lengthscale, we can take 
the limit Ej^ ^ in Eqs. This gives the same 

result as a purely electrostatic calculation [1] . If the dis- 
tance d of the atom from the surface is small, then the 
atom does not feel the curvature of the surface, and one 
expects to get the same energy shift as one would close 
to a plane surface. This is indeed the result we get when 
we take the limit d ~> by using uniform asymptotic 
expansions for the Bessel functions we obtain 

^ 1 „ 1 ^ 1 

=P«^- "**16d^' "^"^IM^- ^^^^ 

2. d<Aji<il 

In this regime the energy shift behaves in exactly the 
same way as in the previous case, because the radius 
of the wire has no influence on retardation, so that the 
relative size of R and Xji does not matter. All that mat- 
ters is that the distance d of the atom from the cylinder 
is still much less than the wavelength A^^ of the rele- 
vant transition in the atom. In mathematical terms, the 



„ 1 ^ ^ ^ 1 

^"^ dHnd' "^'^ d5' ^'"^ dHnd ' 

which is not very helpful numerically though, as logarith- 
mic series converge only very slowly. 



When Xji is smaller than the distance d of the atom to 
the surface of the wire, then the interaction is manifestly 
retarded. As Xji is the smallest of the three lengthscales, 
we first take the limit Xji 0, i.e. Eji 00, in Eqs. 
(P7)) and find that the leading terms in all three integrals 
go as 1/Eji. The remaining integration over k is then 
quite similar to those found in the non-relativistic cal- 
culation in 1] and can be tackled by the same means. 
Scaling fc to a; = kp/m and realizing that the dominant 
contributions to the integrals and sums come from large 
X and large m, one can approximate the Bessel Functions 
by their uniform asymptotic expansions and then gets a 
geometric series, which is easy to sum. In this way one 
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finds the following approximations 
1 



-\Ki{kp)f 



Ax X { yl + x^ H — r= 

1 



K^{kR) 

1 \ A{A^ +AA+1 



(31) 



,3 h{kR) 

2^Ej,pi y Ki{kR) 



P 



dfc k^ 



j dx X 



(A -1)4 
[Ki[kp)Y (32) \, 

A{A^ +AA+1 



(A -1)4 



dx 



Ko(kR) 
3 A{A'^ + iA + l) 



(A-iy 



(33) 



with A{x) given by 



A{x) = 



X exp 




(34) 



These are easy to evaluate numerically and provide a rea- 
sonable approximation to the energy shift in the retarded 
limit, as shown in Fig. ^ In the limit of the distance 
d ~ p — R being much smaller than the radius R of the 
wire, the above approximations yield 



1 1 



(35) 



which agrees with the retarded energy shift of an atom in 
front of a perfectly reflecting plane mirror, as calculated 
by Casimir and Polder [3] . This is what one would expect 
because an atom that is very close to the surface is not 
susceptible to the curvature of the surface. 



\ji-^R<^d 



In this case we again start by taking the limit Eji oo 



in Eqs. ([25)) -(|27 p . which gives a leading order contribu- 
tion proportional to l/Ej^. For distances d much larger 
than the wire radius R the dominant contribution to the 
sum then comes from the summands with the lowest m, 




FIG. 2: The contributions to the energy shift in the retarded 
Umit due to the three components of the atomic dipole, mul- 
tiplied by Ejid'^. Solid lines represent the results of exact 
numerical integration of Eqs. (|25p - (|27|l in the limit Eji — > oo, 
whereas the dashed (red) lines represent the approximations 
(|3ip - (|33p . For large d the asymptotic behaviour is dominated 
by the lowest m terms in the sums, given by (I36|l -- (|38|l and 
shown as dotted (blue) lines. The arrow on the vertical axis 
indicates the exact value in the limit d ^ 0, Eq. (|35p . 



SO that we need consider only those, 



loikR) 



2ttE. 



dkk^4^P^[K,(kp)f 



Ko{kR) 

k I[{kR) 



(36) 



dfc 



p2 K[{kR) 



[Ki{kpt 



1 



i-KEu 



dfcfc fc 



h{kR) 



-2 



ttE 



']l JQ 



dkk 



K[{kR) 

3 /o(fci?) 



Ki{kR) 
K[{kp)] 



[K,{kp)f{ir) 



d,fc3/((fci?) ....... 



K^ikR) 



[Ko{kp)Y 



(38) 



The dotted lines in Fig. [5] show that these are indeed 
good approximations for large d/R. Their leading-order 
behaviour is 



1 



1 



Eji In c? ' 



1 



1 



""^"^ Eu d^' Ej.dHnd ' 



which is in full agreement with the asymptotic results 
by 3, even though those are for a metallic wire char- 
acterized by a plasma frequency. This is because in the 
retarded limit the interaction between the atom and the 
surface depends, to leading order, only on the static po- 
larizability. 

As in the electrostatic case, the contributions due to 
the p and z components of the atomic dipole fall off less 



7 



rapidly than the 4> contribution. We also note that, just 
as in the non-retarded case, the series in powers of 1/ In d 
converge too slowly to be of any practical use, so that 
estimates must be made with Eqs. 

6. i? < Aji < d 

As in the non-retarded cases, the limit of vanishing 
radius {R — s- 0) and the retarded limit {Eji — > oo) com- 
mute, and we recover the results of the previous section, 
Eqs. This is another manifestation of the fact 

that the criterion of whether the interaction is retarded 
depends solely on the distance d between an atom and the 
surface of the wire, and that the relative size of geomet- 
rical features and the wavelength Xji is irrelevant. This 
means in particular that there are no resonance effects 
for Xji coinciding with the wire radius R. 

B. Numerical results 

For intermediate parameter ranges one has to eval- 
uate Eqs. (P5|) - (|?7|) numerically. This is straightfor- 
ward, and one can employ standard software packages 
like Mathematica or Maple. The numerical convergence 
of Eqs. (l25 |) "([27 l) is very good, although more terms are 
needed for small distances d than for large distances. 
Figs. [SHU show the contributions by the p, 0, and z com- 
ponents of the atomic dipole to the energy shift (|15p for 
various values of the typical transition frequency Eji in 
the atom. We give the distance d and the transition 
wavelength 1 / Eji in units of the wire radius R. For plot- 
ting we have factored out of '^p,(t>,z the asymptotic func- 
tional dependence of the shift in front of a plane mirror, 
Eq. (011). 

In Fig. [6] we show how these contributions look when 
we choose the wavelengths 1/Eji of a typical internal 
transition as a lengthscale and plot the contributions to 
the energy shift for various wire radii R. The larger the 
value of R the more terms are required in the numerical 
series. 



IV. ENERGY SHIFT NEAR A PERFECTLY 
REFLECTING SEMI-INFINITE HALFPLANE 

Next we wish to calculate the energy shift of an atom 
in the vicinity of a perfectly reflecting halfplane, as illus- 
trated by Fig. [71 

The procedure of obtaining the normal modes of the 
vector potential is analogous to that described in Section 
mil The scalar solution of the Helmholtz equation ^ in 
the cylindrical coordinates that is best suited to applying 
boundary conditions on the surface of the halfplane is 
given by 
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FIG. 3: The contribution ((25l) to the energy shift ITS)) due to 
the p component of the dipole for various typical transition 
frequencies Eji. The dashed line is this contribution in the 
retarded limit Eji oo. 
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d/R 

FIG. 4: The contribution ((26)) to the energy shift (fTS)) due to 
the (j} component of the dipole for various typical transition 
frequencies Eji. The dashed line is this contribution in the 
retarded limit Eji oo. 



where Jm/2ikp), with m = 0, 1, 2, . . ., are the regular so- 
lutions of Bessel's equation, and the separation constants 
satisfy = k"^ + k^. We must have to > 0, as other- 
wise the solutions are not linearly independent. Note 
that half-integer indices arise because the angle (j) is re- 
stricted to the interval [0, 2n], so that the usual argument 
of single-valuedness of e"""^ cannot be evoked. 

In order to obtain two linearly independent vector so- 
lutions we again apply Eqs. ^ and ©, and impose the 
boundary conditions for a perfectly reflecting halfplane, 
Ep = = Ez and Brj, — for (j) = and cj) = 27r. In this 
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FIG. 5: The contribution ((27)) to the energy shift dTSj) due to 
the z component of the dipole for various typical transition 
frequencies Eji. The dashed hue is this contribution in the 
retarded limit Eji oo. 



way we find for the mode functions 



FW(r) 



\/27r 



^sin(|<^) J^/2(M ep 



Ff(r) 



cos \-T^<^) JL/2{kp) e4 



1 



^/27r 



IK (m \ I , ^ 
— sin y—4>) Jrn/2ikp) ' 



(39) 



(40) 




where the composite index stands for A = {fc, m, k}. For 
m > these mode functions satisfy the normalization 
condition (jlip . but the first polarization has an addi- 
tional mode with m = for which Eq. ([55)1 must be 
multiplied by an additional factor 1 / a/2 for it to be nor- 
malized correctly according to PT|) . 



FIG. 6: The contributions (125}-(127|) to the energy shift ^ 
due to the p, (j), and z components of the dipole for various 
radii R of the wire. 



(41) 




FIG. 7: An atomic dipole in the vicinity of a perfectly refiect- 

incr Gcimi-infJnif Vialfr\lnnci T'Vio nrirmal mrtrloG 



Substituting the mode functions Q-lHIl) into Eq. © 
and renormalizing the energy shift by subtracting the 
free-space contribution in the same way as this was done 
in Eqs. (fTB|) - P^ . we obtain an energy shift of the form 
(fTC)) with 
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dk k dK 



2 oo 



E,,+u; [ykp, ^^^^ 
£'[sin2 {^c^)j'„\„{kp)^j:l{kp) 

( oo 



2 7 V 2 



m=0 



2 /"^ /""^ 
— / dk k dn 
71" Jo Jo 

/ \ 2 oo r 

^ ' ^ m— 1 



, . , COS , 



OO 



(42) 



(43) 
(44) 



r 



where the primes on the sums indicate that the to = 
terms are weighted by an additional factor of 1/2. In 
order to simphfy these expressions, the sums over the 
Bessel functions need to be evaluated. Recently, similar 
summations have been carried out [l6l . [l9| , but the re- 
sults obtained do not include our particular case of sums 
involving Bessel functions of the half-integer order. 

We proceed along the following lines. First, we split 
each sum into two, one over Bessel functions of integer 
orders, and the other over half-integer orders. For the 
first we can apply the standard summation formula (isl . 
9.1.79] 



cos2to0 J^(z) = ■^Jo(2zsin 



(45) 



m=0 



and we choose to represent the right-hand side in terms 
of an integral 9.1.24] 



1 



- Jo(2zsin( 



dt 



sin(2zt sin ( 



1 



(46) 



For the half-integer sum we use a summation formula of 
[20I 5.7.17.(11.)], which in our case gives 



'■1/^'"'^^^ sin(2ztsin0) 



Ecos(2to+1)0J2^^.(z) = 1 / 

m=0 "^1 

(47) 

We note that, if we use the integral representation (j46|l . 
the sums over integer and over half-integer Bessel func- 
tions are very similar; the only difference is the upper 
limit of the t integral in (|46)) and (|47p . As these t inte- 
grals and their derivatives will arise repeatedly, we define 
the following auxiliary functions: 



i^(z,0) 
G{z,4>) 



i/"'"-^^^ sin(2z^sin( 
1 * 



dt 



sin(2zt sin( 



(48) 
(49) 



Further we note that the k integrals in Eqs. (|42p - ((44l) 
suffer from the same convergence problems as already 
discussed in Section IIIII We avoid these by introducing 
polar coordinates with k = wsina and k = cjcosa. At 
the same time we parametrize the denominator arising 
from perturbation theory by 



E. 



dx e-(-^^'+'^)^ 



(50) 



with Ej,+Lu = Eji + Vk^ + k2 > 0. Then Eqs. ((42|) -((44l 
become 



dx e^^^ 



Tr/2 

duj w^e^""^ / da sin a 



X |(Ti((jj/9sina) + (T3 (cjp sin a) cos^ a} , (51) 



dxe"-^^'^ / dujuj^e^'^'' / da sin a 
Jo Jo 

X {(T2(wpsina) cos^ a -|- a4{ujpsma)} , (52) 



da; e"-^^ 



duj I 



IT II 



da sin a 



X {(T5((jj/9sina) sin^ a} 



(53) 



The sums <7i(z) appearing in these expressions can be 
calculated by using Eqs. and standard deriva- 

tive formulae for Bessel functions ^18, 9.1.27]; we obtain 
in terms of ([HI) and ([49]) : 
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(71 (z) 
<T2{z) 



(^;) 
(74 (^) 



m— 1 



sin^ (m0/2) 
cos^ (m0/2) 



87rz2 

OO 

E' 



± 



± 



sin^ (m0/2) 
cos^ {■m(j)/2) 



"V 0=0 



0=0 



(54) 



CTi(z) 

(72 (^;) 



+ ^ [^(^. 0) - G(z, 0)] T [F{z, cj,) + G(z, , 



cos2z, ,, , ^, 

+ ^^(lTcos(/>) , (55) 



OO 



— [F{z, 0) - Fiz, 4>) - Giz, 0) - G(z, 0)] 



(56) 



We now carry out the various integrations in the following 
order. First we evaluate the a integrals, which all give 
Bessel functions Ji or Jq [13, 3. 715(10), (14)]. Next we 
carry out the integrations over to, which involve integrals 
of the type [H 6.611(1.)] 



dz e-^'J^ibz) 



Finally, we calculate the t integrals that came in through 
the auxiliary functions F and G, Eqs. ([35]) and (|49p . 
These are all elementary. At the very end we calculate 
the (j) derivatives of Eq. ([54]) and take the limit (/> — > in 
the appropriate terms. The end results then still contain 
the parameter integral ([50]) over x, which we now scale 
by substituting x = '2p'i]. Then the final results read 



1 



IQttp'^ Jq 



drj e 



[(277^ + 1) sin^ (j>-Tj^ 



(1+772)3/2(^2 



774(1+7^2)3/2 rji (7;2+sin2 0)3 

[(2 + 77^) sin^ (j) + 2sin2 0(37;'' + 677^ + 2) - 77^(37;* + 67?^ + 4)] [> , 



(57) 



1 



167rp3 



drj e 



irf + 67^4 + IO772 + 4 4 



[(1 - 277^) sin^ 4) + 77^] 



cos ( 



(1 + 7y2)5/2(r^2 ^ gjjj2 

POO 



7^4(1+772)5/2 774 (772 + sin^ ,/))3 

[(2 - 2ff - 774) sin4 + 2sin2 0(2 + 27;^ - 67/4 _ 3^6-, ^ ,^2(3^6 ^ g^4 _^ -^q^2 + 4)] | ^ 



1 



167rp3 



d77 e 



-2pEjiri 



9774 + IO772 + 4 4 sin^ 0-772 



(1 + 7y2)5/2(r^2 gjjj2 ^-j3 



774(1 + 7/2)5/2 ^4 (^2_^sin2 0)3 

[(7y2 - 2) sin"* + 2(7/4 _ 4^2 _ 2) sin^ + 77^(9774 + 1077^ + 4)] | 

I 



(58) 



(59) 



Inserted into Eq. ^ the Eqs. give the final 

result for the energy shift of an atom near a perfectly 
reflecting halfplane. Some of the integrations over the 
auxiliary variable 77 could in principle be carried out, but 
those would yield complicated hypergeometric functions. 
Thus it is preferable to have the result in the form of an 
integral over elementary functions. It converges quickly 



and can therefore be very easily evaluated numerically 
by using standard software packages. In addition, we 
shall go on to determine asymptotic expressions in the 
non-retarded and retarded regimes. 
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A. Asymptotic regimes. 

1. Plane-mirror limit 

In the limit of the polar angle (p being very small, the 
atom is very close to the halfplane but far away from the 
edge, so that the energy shift should be the same as for 
an atom in front of a plane, infinitely extended mirror. 
The component of the atomic dipole that is normal to 
the surface should then give the contribution listed in 
Eq. (|29p to the shift, and the parallel components should 
contribute that shown in Eq. ([50]) . As the distance d of 
the atom from the halfplane is psiiKp, we take Eqs. (|57p - 
(|59p and scale rj rj sin ^, so as to get an exponential 
with the same argument as in Eqs. ([M]) and (|5D|) . If we 
subsequently take the limit — > 0, we recover Eqs. 
and ((30)) . as expected. Note, however, that the geometry 
is different from the cylindrical case: the 4> component 
of the atomic dipole is now normal to the surface and its 
contribution to the energy shift is given by (|29p , and 
the p and z components are parallel so that and 
are given by ([501) • 



2. Non-retarded regime 

If pEji <C 1 then the atom is very close to the half- 
plane, compared to the wavelength of a typical internal 
transition. This means that the interaction of the atom 
and the surface is instantaneous, as the atom evolves on 
a much longer timescale. In this case field quantization 
is not necessary, and only Coulomb interactions between 
the atom and the halfplane need to be considered, as was 
done in Ref. [l|, where we derived 



487rp3 
1 



487rp3 
1 



167rp'^ sin 

^ cos (f) 
Sirp^ sin^ 
cos(h 



(tt - 0)(1 + sin^ ( 

IGnp^ sin"^ 
(tt - 0)(1 + cos^ 
IGnp'^ sin^ (j) 

TT — (j) 

167r/9^ sin'^ <j) 



247rp3 167rp3 sin^ 



Taking the limit Eji ^ in Eqs. ([57 j) -(f59 |l we recover 
these results, which is an important consistency check on 
our present calculation. 



distance to the halfplane is p sin (j), and consequently the 
criterion for full retardation is psiiKpEji ^ 1, cf. Fig. [T] 
Taking the limit Eji ^ oo in the integrals ([57| - ((59l) is 
straightforward, since, according to Watson's lemma [2l\ . 
the integral is then dominated by contributions from the 
vicinity of 77 = 0+, so that one just needs to factor out 
the exponential and expand the rest of the integrand in 
the curly brackets in a Taylor series about this point. 
The leading terms of these Taylor expansions turn out to 
be constants with respect to rj in each case. Thus in the 
retarded limit we obtain 



1 



1 



-3 



1 



GAnp'^E. 



sin* (0/2) 
1 

sin^(</>/2) 
1 



sin^(0/2) 
2 



+ 



sin^(0/2) 
2 



sin* (0/2) sin^(0/2) 



(60) 
(61) 
(62) 



which, for the case of isotropic polarizability, is in agree- 
ment with the result of Ref. [1^. In the light of our 
comments above, we emphasize again that these results 
are only valid when the distance of the atom from the 
halfplane exceeds several wavelengths Xji. This means 
that for small angles cj) one needs to revert to the plane- 
mirror limit discussed in Section IIV A II above, because 
in the region < < 7r/2 Eqs. ([501) -(US) apply only if 
sin (j) 3> Xji/ P- However, taking the limit p — s- 00 together 
with 0^0 while keeping psin^ = d fixed is legitimate, 
and reproduces the well-known Casimir-Polder result 
for the retarded interaction between an atom and a plane 
mirror, Eq. ([55)1 . 

Taking the hmit ^ tt in Eqs. ([5(Il) - ([5^ shows that 
for an atomic dipole that is polarized azimuthally the 
interaction vanishes when the atom is located exactly 
above the edge of the halfplane. This conclusion actually 
holds not just in the retarded regime, but generally for 
any distance, as Eq. (|58p also vanishes in the limit (f> — s- 
TT. Purely from symmetry one would expect there to 
be no azimuthal component to the Casimir-Polder force 
directly above the edge, but the fact that there is no 
radially directed force either is surprising. 

Since we have worked in the cylindrical coordinates, 
the direction of the unit vectors Bp and depends on 
the position coordinates p and (p. In this context it is 
curious that, in the retarded limit, all three components 
of the atomic dipole contribute to the energy shift with 
exactly the same angular dependence. 



3. Retarded regime 

In the opposite limit of the atom being far away from 
the halfplane, we need to distinguish whether the atom 
is located beyond the edge of the halfplane or not. If it 
is, i.e. for 7r/2 < < TT the distance of the atom to the 
halfplane is its distance to the edge, namely p, so that 
the condition for the interaction to be fully retarded is 
pEji ^1. If, on the other hand, < < 7r/2 then the 



V. SUMMARY 

We have calculated the energy shift in a neutral atom 
caused by the presence at arbitrary distance of perfectly 
reflecting microstructures of two different geometries. 
For an atom at a distance d = p — i? from the perfectly 
reflecting cylindrical wire of radius R we have found an 
exact expression for the interaction energy, Eq. (|15p with 
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FIG. 8: Direction of the retarded Casimir-Polder force act- 
ing on the atom with isotropic polarizability. Note from Eqn. 
(|58|) that an atom that is polarized azimuthally does not ex- 
perience any force when it is located exactly above the edge 
of the halfplane. 



which was done in detail in Ref. [l|. The case of purely 
retarded interactions, which occur when the distance d 
is much larger than Xji, has been analysed in Sections 
nil A 41 -6. For a small wire radius the three contributions 
to the energy shift are well approximated by Eqs. p6|) - 
(f38l) . and for a large wire radius by Eqs. ([3T|l - (l33|) . 

In the case of an atom close to a perfectly reflecting 
halfplane the exact analytic analysis can be pushed a 
little bit further than in the cylindrical case. We have 
managed to find an exact formula for the energy shift in 
terms of a simple, rapidly converging integral over ele- 
mentary functions, Eqs. (|57ll - ((59|) . so that they are very 
easy to study numerically. Nevertheless, we have also 
derived asymptotic formulae, which agree with previous 
calculations. 

The totality of our results can be used to reliably es- 
timate the energy shift in an atom close to a variety of 
common microstructures that consist of a ledge and pos- 
sibly an electroplated top layer of higher reflectivity. We 
have determined the energy shifts for the complete range 
of distances, which is very important for practical appli- 
cations as in many modern experiments the distance of 
the atom is neither much larger nor much smaller than 
the typical wavelength of an atomic transition, but com- 
mensurate. 



Eqs. ([25|) -(l27 p . As these integrals and sums are in general 
quite complicated, we have analysed various important 
limiting cases. The limit of the distance d being small on 
the scale of the wavelength Xji of a typical atomic tran- 
sition requires only electrostatic forces to be considered, 
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